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Abstract 

Mercury  is  a  50-ns,  6-MV,  360-kA  accelerator  with  a 
magnetically-insulated,  inductive -voltage-adder  (MIVA) 
architecture.  The  machine  was  formerly  known  as 
KALIF-HELIA[1]  at  Forschungszentrum  Karlsruhe  in 
Germany  but  now,  with  some  minor  modifications^], 
will  be  sited  at  NRL.  Mercury  can  be  operated  in  either 
positive  or  negative  polarity[2-4].  Voltage  is  added  in 
vacuum  along  a  magnetically  insulated  transmission  line 
(MITL)  from  six  voltage  adder  cells.  Understanding 
power  flow  and  coupling  to  a  load  in  this  geometry 
requires  the  application  of  MITL  theory[5-8].  Because 
the  electric  field  stresses  on  the  cathode  in  the  MITL 
exceed  the  vacuum  explosive-emission  threshold,  electron 
emission  occurs  and  current  flow  is  divided  between 
current  flowing  in  the  metal  and  in  vacuum  electron  flow. 
This  electron  flow  manifests  itself  as  a  loss  current  until 
the  total  current  is  large  enough  to  magnetically  insulate 
the  emitted  electrons  from  crossing  the  anode-cathode 
(AK)  gap.  Once  insulated,  the  electrons  flow  axially 
toward  the  load  as  illustrated  in  Fig.  1.  In  particular, 
electron  emission  and  flow  along  the  MITL  alters  the 
impedance  along  the  line  and,  thus,  the  power  flow 
coupling  between  the  machine  and  the  load.  The  effective 
impedance  is  best  described  by  the  flow  impedance, 
which  is  a  function  of  both  the  geometry  and  the  voltage. 
When  electrons  are  emitted  from  regions  having  different 
voltages,  such  as  in  the  adders  or  at  different  locations 
along  the  MITL  itself,  layered  flow  occurs,  further 
complicating  the  picture.  Analysis  of  power  flow  in  this 
complex  geometry  is  underway  to  understand  the  past 
performance  of  KALIF-HELIA  and  to  assist  in  optimizing 
the  future  performance  of  Mercury  in  both  polarities  and 
for  various  load  configurations[3,4].  The  goal  of  this 
work  is  to  develop  physics-based  MITL  circuit-element 
models  for  the  NRL  transmission  line  code  BERTHA[9] 
to  properly  treat  power  flow  in  the  vacuum  section  of 
Mercury  while  modeling  the  full  machine. 


I,  MITL  MODEL 

The  fundamental  understanding  of  MITL  flow  is  derived 
from  a  pressure  balance  argument  when  space  charge 
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Figure  1.  Schematic  of  MITL  flow  in  negative  polarity. 


limited  (SCL)  emission  occurs  at  the  cathode[5-8].  This 
assumption  of  SCL  emission  implies  that  the  cathode  is 
turned-on  (i.e.,  the  cathode  is  a  zero  work  function 
emitter)  and  that  there  is  sufficient  space  charge  in  the  gap 
to  drive  the  electric  field  to  zero  at  the  cathode  surface. 
The  model  allows  for  the  possibility  of  additional  electron 
emission  at  any  point  along  the  MITL  if  there  is  not 
enough  space  charge  in  the  gap  and  electron  current  loss 
to  the  anode  if  there  is  too  much  space  charge  in  the  gap. 
It  is  also  assumed  that  there  is  no  ion  emission  from  the 
anode,  the  electrons  are  emitted  from  the  cathode  with 
zero  energy  (and,  hence,  at  zero  pressure),  and  the 
electron  pressure  is  negligible  compared  to  the 
electromagnetic  field  pressure  at  the  anode. 

The  general  form  for  the  flow  impedance  Zf  is 


(1) 


where  Z0  is  the  vacuum  impedance  of  the  line,  V  is  the 
voltage  across  the  line,  Ia  is  the  current  flowing  in  the 
anode,  Ic  is  the  current  flowing  in  the  cathode,  and  Ia  -  It  is 
the  current  flowing  in  the  electron  flow  layer.  This 
definition  for  Zf  is  independent  of  the  distribution  of  space 
charge  in  the  AK  gap.  For  the  case  where  the  space 
charge  density  is  uniform  across  the  flow  layer, 


where  Z0  is  the  vacuum  impedance,  m  and  e  are  the  mass 
and  charge  of  an  electron,  and  c  is  the  speed  of  light.  The 
last  term  is  the  space  charge  correction.  These  two 
equations  can  be  used  to  determine  the  flow 
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Figure  2.  Plot  of  Z0Ia  as  a  function  of  Z0IC  for  various  V. 


characteristics  in  most  situations.  In  the  past  this  result 
has  been  displayed  graphically  by  solving  Eq.  (2)  for  Ia 
and  plotting  Ia  as  a  function  of  Ic  for  various  values  of  V 
as  shown  in  Fig.  2  [5-8].  The  curves  terminate  at  the  left 
when  the  flow  is  saturated  and  the  electron  flow  layer 
extends  all  the  way  to  the  anode.  When  the  flow  is  super 
insulated,  there  is  very  little  flow  current  and  the  flow 
layer  is  very  close  to  the  cathode.  The  minimum  anode 
current  for  a  given  voltage  occurs  for  self-limited  flow, 
which  is  discussed  in  more  detail  below. 

Here,  expressions  for  the  flow  impedance  are  derived 
that  will  be  needed  for  applying  the  results  to  BERTHA. 
When  Z0,  V,  and  Ic  are  known,  Eq.  (2)  is  solved  for  (Ia2  - 
Ic2)12,  which  is  substituted  into  Eq.  (1)  to  obtain 
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Note  that  Zf  =  Zo/2  for  Ic  =  Icins  =  (2mc2V/eZ02) 1/2,  which 
is  the  minimum  cathode  current  required  for  insulation 
(i.e.,  Ic  >  Icms  in  order  for  a  flow  solution  to  exist).  Also 
note  that  Z0/2  <  Zf  <  Z0.  Once  Zf  is  determined,  Ia  is 
found  from  Eq.  (1)  with  Ia  =  (Ic2  +  V2/Zf2)1/2.  Similarly, 
when  Z0,  V,  and  Ia  are  known,  Eqs.  (1)  and  (2)  can  be 


Figure  4.  Plot  of  Zj/Z0  as  a  function  of  Z0Ia  for  various  V. 


manipulated  to  provide 
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where  this  expression  needs  to  be  inverted  to  find 
Zt{Z0,V,Ia).  Here,  Ia  =  Ia0  =  (2V/Z0)(1  +  mc2/2eV)1/2  when 
Zf  =  Zo/2.  Once  Zf  is  determined,  Ic  is  found  from  Eq.  (1) 
with  Ic  =  (Ia2  -  V2/Zf2)1/2.  Plots  of  the  expressions  for 
Zf(Z0,V,Ic)  and  Zf{Z0,V,Ia)  are  shown  in  Figs.  3  and  4, 
respectively,  for  various  values  of  V.  Note  that 
Zt{Z0,V,Ia)  is  double  valued  for  IaSL  <  Ia  <  Ia0,  so  that  care 
must  be  taken  to  determine  which  value  is  appropriate. 
Self-limited  flow  with  Zf  =  ZtSL  occurs  at  the  minimum 
value  of  Ia  =  IaSL  and  describes  the  flow  that  is  established 
in  a  long  MITL  or  an  MITL  coupled  to  a  high  impedance 
load  (e.g.,  for  ZL  >  ZfSL). 

Analytic  expressions  for  the  flow  impedance  and 
currents  in  self-limited  flow  have  been  derived  for 
application  in  BERTHA  and  are  given  by 

=  ZoIsl(V )  >  where  (5) 


SL 

^  f 


fSL(V)  = 


Vn  ^ 

f vn  } 

4V 

± 

4V 

Z-i 

VSL-, 

V  J 

V  t 

~\l  /  2 


Vo. 

V 


-1 


(6) 


r 


if  (V  )  = 


If(V)= 


Z0fSL(V) 


1 


,  1/2 


2Z 


o  J 


-1  + 


\  If st,  (V ) ' 

20V 


,  and  (7) 


(8) 


o 


+ 


,  12V  36V-  32V3 

1  + - + - ^  + - r- 


“1 1/2  \l/2 


Here,  V0  =  mc2/2e  =  0.255  MV,  and  in  Eq.  (6)  the  upper 
(+)  sign  applies  when  V0/V  >  1  and  the  lower  (-)  sign 
applies  when  V0/V  <  1.  Note  that  Zt/Z0,  Z0Ia,  and  Z0IC  are 
functions  only  of  V 
are  shown  in  Fig.  5. 
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Figure  5.  Plots  of  ZfSL/Z0,  Z0IaSL,  and  Z0ICSL  as  functions 
of  V. 


II.  BENCHMARKING  MODEL  RESULTS 
WITH  PIC  SIMULATIONS 

Numerical  simulations  of  simple  MITL  flow  problems 
have  been  run  to  benchmark  the  model.  The 
electromagnetic  hybrid/PIC  code  LSP  [10]  was  used  to 
ran  these  simulations.  A  long  (compared  with  the 
distance  required  to  establish  equilibrium  flow)  MITL  that 
is  terminated  with  an  electron  beam  diode  load  is  driven 
in  negative  polarity  (i.e.,  the  center  conductor  is  the 
cathode)  to  equilibrium  at  a  voltage  of  1  MV.  The 
simulation  region  is  10  cm  long  with  a  cathode  radius  of 
0.5  cm  and  an  anode  radius  of  1  cm,  which  represents  a 
vacuum  line  impedance  of  41.6  Q.  The  model  predicts  a 
self-limited  flow  impedance  of  Zf  =  32  O  for  this 
geometry  at  1  MV  (see  Eqs.  (5)  and  (6),  or  Fig.  5).  Three 
cases  were  run  with  axial  diode  gaps  of  D  =  0.25  cm,  0.53 
cm,  and  1  cm.  Particle  plots  of  these  three  cases  are 
shown  in  Fig.  6.  At  D  =  0.25  cm,  the  diode  impedance  is 
about  15  Q,  well  below  ZfSL.  In  this  case,  the  equilibrium 
flow  is  on  the  super- insulated  branch  of  Fig.  4  with  the 
flow  layer  held  close  to  the  cathode  as  seen  in  Fig.  6a. 
The  currents  measured  in  the  simulation  at  z  =  5  cm  are  Ia 


z(cm) 


Figure  6.  Plots  of  MITL  electron  flow  driven  at  1  MV 
with  (a)  D  =  0.25  cm,  (b)  D  =  0.53  cm,  and  (c)  D  =  1 


=  48.2  kA  and  Ic  =  40.4  kA,  which  corresponds  to  a  flow 
impedance  of  37.7  Q  as  calculated  from  Eq.  (1).  Based 
on  Ia  =  48.2  kA,  Eq.  (4)  predicts  a  flow  impedance  of  37.4 
Q  and  a  cathode  current  of  40.3  kA,  which  are  in 
reasonable  agreement  with  the  simulation  results.  This 
implies  that  the  assumption  of  uniform  charge  density  in 
the  flow  layer,  which  was  used  to  derive  Eq.  (2)  (and  the 
subsequent  equations),  is  reasonable.  The  charge 
densities  measured  in  the  simulations  are  compared  with 
the  assumed  charge  densities  in  Fig.  7.  Although  the 
charge  density  profiles  found  in  the  PIC  simulations  are 
peaked  near  the  cathode  and  contain  somewhat  more 
space  charge  compared  with  the  square  profiles  assumed 
in  the  model,  the  widths  of  the  charge  layers  are  in 
reasonable  agreement. 

For  D  =  0.53  cm,  the  load  impedance  is  comparable  to 
the  self-limited  flow  impedance.  In  this  simulation  shown 
in  Fig.  6b,  Ia  =  44.1  kA,  Ic  =  29.7  kA,  and  Zf  =  30.6  Q  at  z 
=  5  cm,  which  are  approaching  to  the  self-limited  flow 
values  of  IaSL  =  42.5  kA,  ICSL  =  28.9  kA,  and  ZfSL  =  32  Q. 
predicted  by  the  model.  For  D  =  1  cm  as  shown  in  Fig  6c, 
the  load  impedance  is  significantly  larger  than  the  ZfSL 
and  the  MITL  is  expected  to  run  self-limited.  The 
simulation  values  of  Ia  =  43.1  kA,  Ic  =  28.5  kA,  and  Zf  = 
30.6  Q  are  in  reasonable  agreement  with  the  predicted 
self-limited  flow  values. 


Figure  7.  nc(r)  from  PIC  simulations  with  (a)  D  =  0.25 
cm,  (b)  D  =  0.53  cm,  and  (c)  D  =  1  cm. 
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III.  INCORPORATING  THE  MITL 
MODEL  INTO  BERTHA 

The  TRIFL  code  was  developed  previously  to  model 
the  MITL  flow  of  pulsed  power  vacuum  transmission 
systems[  11].  The  purpose  of  the  work  presented  here  is  to 
incorporate  the  MITL  model  into  a  transmission  line  code 
in  order  to  develop  the  capability  to  model  a  complete 
system  from  the  primary  energy  store  to  the  load 
including  an  MITL  section.  The  transmission  line  code 
BERTHA  will  be  used  for  this  purpose[9].  An  algorithm 
for  incorporating  the  MITL  model  into  BERTHA  needs  to 
be  flexible  enough  to  handle  a  number  of  important 
processes.  The  MITL  is  modeled  as  a  series  of  short,  one- 
time-step-long,  variable  impedance  transmission  line 
elements  each  with  a  variable  shunt  resistance.  The 
impedance  of  each  element  is  set  at  the  instantaneous 
local  flow  impedance  and  the  shunt  resistance  is  adjusted 
to  account  for  losses  in  the  insulation  front,  at  impedance 
transitions,  and  in  the  load-coupling  region.  Where  there 
are  no  losses,  the  shunt  resistance  is  set  to  a  large  value. 

The  algorithm  must  operate  in  either  positive  or 
negative  polarity  and  account  for  voltage  reversal,  while 
handling  cathode  turn-on  and  emission  in  the  insulation 
front,  current  loss  to  the  anode  in  overmatched  impedance 
transitions,  addition  to  the  flow  current  in  undermatched 
impedance  transitions,  voltage  addition  at  adder  junctions 
with  and  without  flow  in  the  adders,  self-limited  flow 
when  appropriate,  re-trapping,  and  flow  coupling  to  a 
load.  At  each  time  step,  the  voltage  waves  on  the 
transmission  line  are  advanced  and  then  the  impedance 
and  shunt  resistance  of  each  element  are  adjusted  to 
account  for  the  change  in  flow  conditions. 

Once  turn-on  occurs,  electron  emission,  limited  to  the 
SCL  current  density,  launches  flow  into  the  MITL 
section.  The  flow  impedance  is  calculated  for  each  MITL 
element  starting  from  the  upstream  end  of  the  MITL 
working  downstream  to  the  load.  When  a  downstream 
MITL  element  is  overmatched  to  its  upstream  neighbor, 
current  is  lost  to  the  anode  but  the  cathode  current  is 
preserved  across  the  transition.  In  this  case,  V  and  the 
downstream  Ic  are  known  and  Eq.  (3)  (or  Fig.  3)  is  used  to 
find  the  downstream  Zf.  When  a  downstream  MITL 
element  is  undermatched  to  its  upstream  neighbor, 
additional  current  is  launched  into  the  flow  but  the  anode 
current  is  preserved  across  the  transition.  In  this  case,  V 
and  the  downstream  Ia  are  known  and  Eq.  (4)  (or  Fig.  4)  is 
used  to  find  the  downstream  Zf.  Matching  conditions  at 
transitions  where  adders  connect  to  the  MITL  have  also 
been  derived  in  order  to  apply  Eq.  (3)  or  (4)  at  these 
junctions.  At  the  load,  the  flow  current  combines  with  the 
load  current  or  is  lost  to  the  anode  depending  on  the  load 
properties.  For  low  impedance  loads  the  flow  in  the 
MITL  is  determined  by  the  load  impedance,  while  for 
high  impedance  loads  self-limited  flow  is  established. 

This  algorithm  for  incorporating  the  MITL  model  into 
BERTHA  has  been  developed  and  is  being  implemented. 


IV.  SUMMARY 

A  model  has  been  developed  to  describe  MITL  flow  in 
a  vacuum  transmission  line  section  of  pulsed  power 
generators  for  implementation  in  the  transmission  line 
code  BERTHA.  The  immediate  application  will  be 
analysis  of  power  flow  in  the  new  Mercury  accelerator  at 
NRL,  which  is  based  on  MIVA  technology.  The  model 
has  been  successfully  benchmarked  against  basic  PIC 
simulations  and  the  algorithm  for  applying  the  model  is 
now  being  incorporated  into  BERTHA. 
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